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ABSTRACT 

Stellar multiplicity lies at the heart of many problems in modern astrophysics, including the physics 
of star formation, the observational properties of unresolved stellar populations, and the rates of 
interacting binaries such as cataclysmic variables, X-ray binaries, and Type la supernovae. However, 
little is known about the stellar multiplicity of field stars in the Milky Way, in particular about the 
differences in the multiplicity characteristics between metal-rich disk stars and metal-poor halo stars. 

In this study we perform a statistical analysis of ^14,000 F-type dwarf stars in the Milky Way through 
time-resolved spectroscopy with the sub-exposures archived in the Sloan Digital Sky Survey. We obtain 
absolute radial velocity measurements through template cross-correlation of individual sub-exposures 
with temporal baselines varying from minutes to years. These sparsely sampled radial velocity curves 
are analyzed using Markov chain Monte Carlo techniques to constrain the very short-period binary 
fraction for field F-type stars in the Milky Way. Metal-rich disk stars were found to be 30% more 
likely to have companions with periods shorter than 12 days than metal-poor halo stars. 

Subject headings: binaries: close — binaries: spectroscopic — Galaxy: stellar content — stars: statis¬ 
tics — surveys 


1. INTRODUCTION 

Stellar multiplicity plays a crucial role in many fields 
of astronomy. Star formation and evolution, Galactic 
chemical evolution, nuclear astrophysics, and cosmology 
are all influenced by our understanding of the multiplic¬ 
ity properties of an underlying stellar population. Binary 
interactions lead to phenomena as diverse as cataclysmic 
variables, classical novae, X-ray binaries, gamma-ray 
bursts, and Type la supernovae. Stellar interactions are 
also the cause of the anomalous surface abundances mea¬ 
sured in Ba stars, CH stars, and the majority of c arbon- 
enhanced metal-poor stars dLucatello et al.ll2005D . The 
rates of these phenomena depend on the multiplicity 
properties such as the fraction of stars with compan¬ 
ions and the distributions of separations and mass ra¬ 
tios. How these properties are in turn affected by vari- 
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ables such as stellar age, metallicity, a nd dynamic al envi¬ 
ronment remains poorly understood. iMoe fe Di Stefana 
(1201311 find no significant trends with metallicity for O- 
and B-stars, but more work is needed for lower-mass 
stars. 

The recent review by iDuchene &: Krausi (12013T) sum¬ 
marizes the state of the art in multiplicity studies. 
The fraction of systems with companion s is known 
to be a strong func t ion of stellar ma ss (lLadal 12003: 
iRaghavan et al. 120101: iClark et~all I2012T) . and there are 
hints that lower mass systems have smaller separations 
(IDuquennov fe Mavodll99lt IAllenll2007t IRaghavan et al.1 

1201011 . Studies of the Solar neighborhood also indicate 
that lower metalli city stars are more likely to have stel¬ 
lar companions (Rag havan et al.ll2010ll . 

These results are based on heterogeneous samples of a 
few hundred stars at most, often dominated by wide sys¬ 
tems which will never become interacting binaries. The 
spectroscopic surveys that reach small periods are labor 
intensive because large numbers of radial velocities (RVs) 
are required to find the orbital solution of each target. 
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This leads to small sample sizes, which have only in¬ 
creased modestly in the pas t two deca des, from 167 in 
iDuauennov fe Mavod (|1991fl to 454 in iRaghavan et al.1 
(|2010f) . The drive to collect complete samples has limited 
previous spectroscopic studies to the Solar neighborhood 
or specific stellar clusters, but neither of these strategies 
can probe the full range of metallicities and ages span¬ 
ning the field stars of the Milky Way (MW) disk and halo 
components. These limits bias the interpretation of data 
against the global properties of, and variation within, the 
MW field. Thus, we are motivated to take a statistical 
approach with a sample of stars located throughout the 
field in order to investigate their multiplicity properties 
with respect to age, [Fe/H], and component membership. 

With the advent of multip lexed spectroscop ic surveys 
like SDSS (|York et al.l 12000( 1 and LAMOST (jCui et, al.1 
[Ml), we can use multiple RV measurements of thou¬ 
sands of stars to study the properties of stellar multi¬ 
plicity that are more r epresentative o f the entire Galaxy. 
SDSS Data Release 8 (lAihara et al.l [20TTh contains over 
1.8 million optical spectra from the original SDSS spec¬ 
trographs including over 600,000 stellar spectra. In 
this work we employ a lesser known SDSS feature, the 
time-resolved dimension. To facilitate cosmic ray re¬ 
moval, spectra were constructed through co-addition 
of several individual sub-exposures, typically 15 min¬ 
utes in duration. Although under-utilized, the bene¬ 
fit of t he sub-exposu r e dom ain is recognized in work s 
such as lBadenes et, al.l (|2009f) and lBickerton et all (12012th 
Portions of the sky were also re-observed for calibra¬ 
tion and scientific purposes. These additional pointings, 
combined with the sub-exposures, yield a time dimen¬ 
sion where single stars have exposure coverage ranging 
from 3 sub-exposures up to over 40 sub-exposures, and 
time gaps from hours to nearly a decade. The tech¬ 
niques employed herein follo w the time-resolved work by 
iBadenes &; Maod (120121) and lMaoz et al.l (1201211 . 

2. MEASUREMENTS 

2.1. SDSS Observations and Sample Selection 

F-type dwarfs are chosen for our sample because of the 
large number of stars targeted by SDSS with repeat ob¬ 
servations, and their relatively mild variability and ac¬ 
tivity. Additionally, F-stars have main sequence (MS) 
lifetimes greater than 5 Gyr, allowing us to select MS 
stars from both the younger disk and ol der halo. The 
Sloan Stellar Parameter Pipeline !SSPP: lLee et al.l[2008f) 
was developed to determine parameters for stellar spec¬ 
tra in the SDSS archive, including metallicity [Fe/H], ef¬ 
fective temperature T e ff, and surface gravity log g. Sam¬ 
ple selection began with i dentify ing science primary ob¬ 
jects from SEGUE-1 (lYannv et al. l200Dh and SEGUE-2 
(Rockosi et al., in prep.) in the SSPP that were clas¬ 
sified as an F-type star by the “Hammer” classification 
code (jCovev et al.ll2007T ). To minimize the effects of stel¬ 
lar evolution on multiplicity, we selected only dwarf stars 
(log g > 3.75). Stars with multiple fiber pluggings were 
identified astrometrically and joined with the appropri¬ 
ate science primary fibers. 

After measuring stellar RVs ISection 12.21) . systematics 
were revealed in the SDSS sub-exposure spectra. These 
correlations appear as similar shifts in RVs for many 
fibers located on the same plate, typically affecting neigh¬ 


boring fibers on the CCD. After plate-wide comparisons 
of F-stars, RV correlations were corrected where possi¬ 
ble. Corrections applied to the 10 4 RVs are as large as 
17 km s' 1 with a standard deviation of 2.2 km s _ . Not 
all correlations could be identified automatically because 
of multiple groups of correlated shifts, opposite in direc¬ 
tion, on some plates. Visual inspection of plates contain¬ 
ing numerous false binary detections lead to the removal 
of 25 plates including 1155 stars. We urge individuals us¬ 
ing sub-exposure spectroscopy in SDSS to consider these 
systematic shifts in the wavelength solutions. 

Quality control consisted of the removal of: stars with¬ 
out valid parameters in SSPP, fibers located on ‘bad’ 
plates, sub-exposures with a median pixel signal-to-noise 
ratio (SNR) less than 20 or with fewer than 3000 un¬ 
flagged pixels, stars with time lags AT < 1800 s, stars 
with less than three clean sub-exposures, and corrupt or 
misclassified spectra (from visual inspection of stars with 
the largest RV variation or non-characteristic T e s). The 
final sample consists of 14,302 stars (16,894 fibers) with 
as many as 47 sub-exposures, spanning up to nine years 
of observations (Figure [H). 

Our cleaned sample is characterized by metallicities 
ranging from —3.41 < [Fe/H] < +0.52. To aid com¬ 
parison in our analysis, the final sample was sub-divided 
into three groups of equal size by cuts in metallicity at 
[Fe/H] = -1.43 and [Fe/H] = -0.66 (Figure [[]). The 
majority of the stars have three or four sub-exposures 
(median = 4), typically taken about 15 minutes apart. 
The median time lag for a star is 2 hours, however more 
than three years between observations can be seen in 
more than 250 stars (Figure Q]). SNRs for sub-exposures 
lie in the range 20 < SNR < 84 with a median value of 
32. 

2.2. Radial Velocities 

RV measurements were attained through cross¬ 
correlation of sub-exposures with a master template con¬ 
structed from 7207 sample-star, co-added spectra where 
the co-added SNR > 50. The spectra were de-shifted 
using the redshift value assigned to the co-adds by the 
SDSS pipeline, continuum-normalized, and averaged to¬ 
gether. 

Sub-exposures were independently prepared and cross- 
correlated with the template. Spectra were continuum- 
normalized by dividing the spectrum with a highly 
smoothed version of itself using a FFT smoothing al¬ 
gorithm, and then cross-correlated with the template at 
various integer pixel lags. Each spectrum had a cross¬ 
correlation function (CCF) that was fit with a smooth 
spline interpolation. With spectral resolution of R ~ 
2000, the peak lag in pixels translates to the spectrum’s 
redshift at 70 km s _1 pixeU 1 . The mean and standard 
deviation of RVs for individual stars are shown in the Fig¬ 
ure [2] distributions. The velocity dispersion of the mean 
RVs decreases with increasing [Fe/H], indicating that our 
[Fe/H]-groups sample both the disk and halo components 
of the MW. The standard deviation of RVs within indi¬ 
vidual stars is larger for the metal-poor group; however, 
empirically estimated uncertainties also show larger mea¬ 
surement errors for metal-poor stars. This underscores 
the importance of the use of proper error analysis in a 
method such as ours. 
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Fig. 1.— Left: Metallicity distribution for 14,302 F-dwarfs. Right: Distribution of maximum time lag between the first and last exposure 
of a star. 



Fig. 2.— Mean (left) and standard deviation (right) of radial velocities within a star. Variations in the standard deviation of velocities 
are affected, in part, by the larger measurement uncertainties for metal-poorer stars. 


2.3. Uncertainties 


It is well known that uncertainties in CCF peaks must 
be estimated empiricall y or t hrough some Monte Carlo 
method (e.g., iPeterson et al.lfl998l) . For this work we 
determined RV uncertainties empirically by quantifying 
the spread in measurements for spectra of similar qual¬ 
ity. The median absolute deviation (MAD) is a ro¬ 
bust measure of the variability of a sample and is re¬ 
lated to the standard deviation by cr = 1. 482 6MAD, 
where MAD = median(|RVi — median(RV)|) dLevs et al'l 
120131 ). All measurements were de-shifted into the rest 
frame using the SDSS estimates of the co-add redshift, 
and placed into bins of similar metallicity ([Fe/H] ±0.25) 
and signal-to-noise (SNR ± 2.5). Initial tests showed 
no correlations between measurement spreads and either 
log g or T e ff. Estimates for the uncertainty of RV mea¬ 
surements within a bin were calculated using MAD val¬ 
ues. Here, it is assumed that the majority of stars do not 
have detectable variability over the observed time base¬ 
line, and that effects from intrinsic variations in RV are 
minimized by adopting median values. After performing 
this process for all bins, a functional form for assign¬ 
ing RV measurement uncertainties <r R v was fit with an 
inverse proportionality to SNR, and with a linear cor¬ 
rection in [Fe/H]. The measurement uncertainty as a 
function of [Fe/H] and SNR is, in km s~ , 


a RV ([Fe/H], SNR) 


(-26.51 [Fe/H] ±50.52) 
SNR 


1.23. 

(1) 


Uncertainties are sub-pixel, falling below the spectral res¬ 
olution of 70 km s -1 pixeU 1 . For exposures with SNR < 
25, uncertainties range from 3.0 to 8.0 km s~ , with a 
median value of 5.0 km s . Exposures with SNR > 40 
have uncertainties in the range 1.9 to 4.4 km s _ , with a 
median value of 2.7 km s -1 . 

3. MULTIPLICITY 

The probability of a star having a companion was 
determined through model comparison using a trans- 
dimensional, hierarchical, Markov chain Monte Carlo 
(MCMC) method. Two models were compared: a single¬ 
star model A4 S , and a binary-star model A 4b- The hyper¬ 
parameter A, indexes the model choice at each step in the 
MCMC chain. We evaluated the hierarchical model using 
the Python pac kag e emcee, a MCMC ensemble sampler 
([Foreman-Mackev et al.ll2013i ). 

The single-star model A4 S , fits a star with non-varying 
RVs, parameterized by a systemic velocity Vq. Because 
intra-plate systematics are known to exist, it is reason¬ 
able to assume inter-plate systematics exist as well. In 
light of this, (P — 1) additional parameters pSi, were in¬ 
cluded for each star, where P is the number of plate-MJD 
pluggings composing the star. These plate-shift param¬ 
eters allow all RVs from plate i to shift by some amount 
pSi , relative to the first plate Pq. For the majority of 
stars P = 1, no plate-shift parameters are necessary, and 
A4 S is a 1-parameter model. 

In the binary star model A 4b, the sparsely sampled 
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RVs are fit by a sinusoid defined by four-parameters: the 
log of the semi-amplitude log A, the log of the period 
logP, the phase f>, and the systemic velocity Vq. We 
assume circular orbits (eccentricity, e = 0 ), which is a 
safe assumption fo r tidally circularized, s hort-period or¬ 
bits (P < 12 days: iRaghavan et al.l 12010 ). where we are 
most sensitive. A small number of the binaries found 
in this study may have longer periods and could have 
non-zero eccentricities, but this does not affect our re¬ 
sults. Plate-shift parameters were also adopted in Mb 
wherever P > 1. 

Uninformative priors were used in the MCMC. The 
model index A, has a flat prior from 0 to 1, where 
A < 0.5 denotes M s and A > 0.5 denotes Mb- The 
semi-amplitude prior is log-uniform from 3 km s _ , com¬ 
parable to the measurement uncertainties where M s and 
Mb become degenerate, to 250 km s~ , greater than the 
largest RV differences in the sample. The prior on the 
period is uniform in the range 4 < logP (s) < 7. The 
lower limit logP (s) = 4.0 is equal to the orbital pe¬ 
riod at which stellar contact is certain for low-mass com¬ 
panions. Above logP (s) = 7.0, RV amplitudes in bi¬ 
nary systems are comparable to the measurement uncer¬ 
tainties. Combined with the sparsity of the RV data, 
systems with periods longer than logP (s) = 7.0 are 
outside our range of sensitivity. Priors are also uni¬ 
form for the phase (0 < <f> < 2ir) and systemic velocity 
(—600 < Vo (km s~ ) < 600). Markov chains were run 
independently on every star with an ensemble of 200 par¬ 
allel chain “walkers” for a total of 2.4 x 10 6 samples, then 
burned and thinned to 6 x 10 5 independent samples of 
the posterior. 

Evidence for detection of a companion star is reflected 
by the relative probabilities of A. We define the proba¬ 
bility for the binary model, 77 as the fraction of samples 
in the marginalized posterior having A = Mb- We note 
that the value of 77 is dependent on the choice of priors, 
and is sensitive to the treatment of the SDSS systemat- 
ics. Moreover, a degeneracy arises as the RV curve of 
a long-period, low-amplitude system becomes indistin¬ 
guishable from a single-star system. With this mind, we 
stress that values for 77 are not absolute probabilities of a 
system having a companion, but reflect the ability of the 
data to rule out models under the given prior. However, 
the [Fe/H]-groups can be compared, relatively, by con¬ 
sidering the fraction of systems where 77 is large and M s 
is strongly disfavored. The results are shown in Figure [3] 

We also i n vest igated the e/i parameter proposed by 
iGeller et al.1 (1200811 as a metric for identifying the stars 
with large RV variations. We find that the e/i parameter 
singles out many of the same stars as our more sophis¬ 
ticated MCMC-based inference. Our method not only 
takes into account deviations in RV from the mean, but 
also how well the data fit the expected periodicity of a 
binary system. 

Analysis of the posterior, and visual inspections of the 
binary model fits, show that 681 stars with 77 > 0.65 
are true spectroscopic binaries, though given the spar¬ 
sity of the RV curve sampling, there are sometimes large 
uncertainties in the fitted values for specific model pa¬ 
rameters. Another natural break point is 77 > 0.95; these 
are 209 stars for which the determination and analysis of 
accurate individual model parameters should be possible 


(and will be characterized in future work). An intermedi¬ 
ate cut at 77 > 0.80 is a compromise between these limits, 
yielding a larger sample of stars (406) with modest model 
constraints. The values of the binary fractions that we 
derive below are insensitive, within the uncertainties, to 
the exact choice of cut in 77 . This implies that the RV 
variations for our binary detections are sufficiently above 
the measurement uncertainties, and that the binary frac¬ 
tions reported are not biased due to differences in SNR 
or absorption features. 

Figure [4] shows the logP posteriors for each [Fe/H]- 
group, marginalized over all binary systems (77 > 0.80). 
The posterior distributions of log P for many of these 
stars are complex: many are multimodal, affected by 
aliasing or other issues related to the sparse, biased time 
sampling. One such effect is the increase in probabil¬ 
ity at logP = 4. Here the metal-rich and -intermediate 
groups contain more stars than the metal-poor group 
with At ~ 10 4 s. Systems with periods as short as this 
are extremely rare dDrake et al.ll2014f) . and our increased 
probability in this area may be due to overfitting. Addi¬ 
tionally, the gap at At = 10 4 6 s = 12 h (Figure [TJ may 
affect the estimate of a period. We defer a more sophisti¬ 
cated analysis to a future paper, but these effects should 
not alter the ability to rule out a single-star model. For 
now, Figure 0] illustrates that we are mainly sensitive to 
periods in the range 4 < log P (s) < 6 , or less than about 
12 days. We emphasize that a more detailed analysis will 
be necessary to estimate the true underlying logP dis¬ 
tribution in our sample. 



Fig. 3. — Distribution of 77, the fraction of posterior samples using 
the binary model, for stars. 

4. DISCUSSION 

In Figure [5] we show fb , the measured lower bound 
for the fraction of stars with short-period companions 
(P < 12 days) for each metallicity group, normalized to 
the metal-rich binary fraction, fb is a lower limit because 
of non-detections as a result of sparsely sampled RVs 
and high orbital inclinations, resulting in low amplitudes. 
We see agreement in fb measured for all three choices in 
77 cutoff (0.65, 0.80, 0.95). With a cutoff of 77 = 0.80, 
values of fb for the metal-poor, -intermediate, and -rich 
groups respectively are: 2.5% ± 0.2%, 2.8% ± 0.2%, and 
3.2% ± 0.3%. Since the observational biases that affect 
binary detection are mostly due to the sparsity of the RV 
coverage, which isn’t metallicity-dependent, we conclude 
that the field F-type MS stars in our metal-rich sample 
are, at a 2-sigma level, 30% more likely than those in our 
metal-poor sample to have close binary companions. 
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Fig. 4. — Averaged probability distributions of log P for all binary 
detections (r] > 0.80). These do not reflect actual distributions of 
periods, and should only be used as a guide to probe the region of 
MCMC sensitivity. The shaded region indicates where Roche lobe 
overflow and contact becomes relevant. The dashed line marks the 
circularization limit at a period of 12 days. 

Our metal-rich and metal-poor samples mostly trace 
the MW disk and halo. Differences in the fraction of 
short-period systems can stem from differences in the 
star-formation process, dynamical interactions after star 
formation, or some combination of the two. 

Three-dimen si onal hydrodynamic models from 
iMachida et al.l (|2009f ) actually suggest a higher fre¬ 
quency of binaries formed through cloud fragmentation 
for metal-poor clusters, due to the decreased require¬ 
ment of a cloud’s initial rotation energy to fragment. 
Moreover, their models yield systems with shorter initial 
separations at lower metallicities. The increased fb 
observed for metal-rich stars in this work can more likely 
be explained by dynamical processes than by formation 
processes. 

The observed differences in fb could be explained if 
the clusters that yielded halo field stars had larger stel¬ 
lar densities an d /or g as densities than those of the disk. 
iKorntreff et al.l (|2012l) explore the effects of gas-induced 
orbital decay on period distributions in clusters. They 
note that an increased density of gas in a newly formed 
cluster will lead to a larger number o f short-period sys - 
tem mergers shortly after formation. iParker et al . ; 2009 ; 
describe how clusters with higher stellar densities destroy 
wide binaries through dynamical interactions. An in¬ 
crease in the destruction of high-mass, wide-binary sys¬ 
tems leads to the ejection of former F-star secondaries 
into the field. These orphaned, single-star systems would 
increase the total number of F-star systems in the halo 
field, effectively decreasing the short-period binary frac¬ 
tion measured. Observational evidence of these denser 
cluster environments is needed to support these argu¬ 
ments for a lower fb in the halo. 

Additionally, some close binaries may also transfer 
mass and covert themselves into blue stragglers dLu et al.l 
I201(T ) . Evidence for a n abundance of blu e stragglers in 
the halo has been seen dYannv et al.ll2000h . and may con- 
tribu te to the lower fb observe d in the metal-poor group. 
Also. lDuchene &; Krausl (l2013i) show a decrease in fb with 
age for Solar-type stars, although this result is based on 
visual binaries with wider periods, and is poorly con¬ 
strained due to limited sample sizes._ 

We note that th e recent results of I Gao et al.l (|2014T) and 
lYuan et all (120141 ) . using data from SDSS, show a larger 
binary fraction for metal-poor than metal-rich FGK stars 
in the field. In addition to probing longer periods, the 


former work does not make use of sub-exposure infor¬ 
mation (using only two RV epochs per star) and relies 
on the correctness of model values for the period distri¬ 
bution, mass ratio distribution, and initial mass func¬ 
tion. The latter work, which uses photometric color 
deviations to infer companions, shows a modest metal- 
licity dependence on total binary fraction. Since their 
method is not sensitive to period, the binary fractions 
they report are strongly dominated by more common, 
wider-period systems near the peak of a log-normal pe¬ 
riod distribution (logP (s) = 10 for nearby, Solar-like 
stars: [Raghavan et al. 2010). It is clear that conclusions 
about binary fraction depend on a number of factors, 
especially the range of periods to which the search is 
sensitive and assumptions made about the overall period 
distribution. 

Our MCMC analysis yields posterior probabilities in 
parameter space, allowing for a more detailed study of 
binary properties (e.g., period and separation distribu¬ 
tions) , which will be presented in future work. The tech¬ 
niques in this work have direct applications for current 
and future multiplexed spectroscopic surveys. 
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[Fe/H] 

Fig. 5. — Short-period binary fraction limits, relative to the metal-rich group. Binary companion detections are defined by a cut in 77, 
the fraction of posterior samples using the binary model. Group median values of [Fe/H] are used. 
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